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A partial wave analysis of the system produced in the charge exchange reaction: 

-K~p —^ TT^TT^n at an incident momentum of 18.3 GeV/c is presented as a function of Tr'V 
invariant mass, m^o^o, and momentum transfer squared, from the incident to the 
outgoing tt^tt" system. For small values of the 5- wave intensity shows a broad enhance- 
ment at low m^o^o with a sharp dip in the vicinity of the /o (980). A dip is also observed in 
the vicinity of the fo (1500). There is rapid variation of the 5* to Do relative phase difference 
in these mass regions. For large values of \t\, the fo (980) appears as a bump. The /2 (1270), 
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observed in the D- waves, is produced dominantly by yr-exchange at low values of \t\ and 
02 -exchange at higher values of |t|. 

13.25.Jx, 14.40. Cs 



I. INTRODUCTION 

This paper reports on a high-statistics partial wave analysis (PWA) of the vr^Tr" system produced in 
the charge exchange reaction: Tr~p —^ 71*^71*^71 at an incident momentum of 18.3 GeV/c using data taken 
by experiment E852 at Brookhaven National Lab (BNL). The PWA was performed over the mT^OT^o mass 
range from near threshold {2m^o) to as high as 2.2 GeV/c^ in 0.04 GeV/c^ mass bins and in various bins in 

momentum-transfer-squared t — ~ Ptt-kI^ ~ \Pn ^ Pp\^ ■ 

The mass and \t\ dependence of tttt production in n^- induced reactions with one pion exchange (OPE) 
provides information on the process tttt — > tttt, involving the scattering of the lightest hadrons The 
extraction of vrvr — > tttt amplitudes is, however, complicated by the presence of production mechanisms other 
than OPE ^,||. The \t\ and rriTTTr-dependence of the partial wave amplitudes and their relative phases, the 
focus of this paper, provide information on these mechanisms and the necessary input for future tttt scattering 
studies. 

The study of the tttt system also bears on current issues in the spectroscopy of conventional qq mesons 
and non-gg mesons such as glueballs or mesonic molecules. In particular, the isoscalar scalar and tensor 
sectors have more states than can be accommodated within the conventional qq model. A recent review 
of light meson spectroscopy Q includes a summary of the current experimental situation in these sectors. 
Non-gg candidates include the poorly understood /o(980) and the glueball candidates /o(1500) and /./(ITIO), 
all of which couple to the 7r7r system [Q. Information about the masses, widths, and decay modes of these 
states, along with knowledge of their production mechanisms, as revealed by their \t\ dependences, will help 
in unraveling their substructure |ll|-|lq| . A complete understanding of these states requires corresponding 
information from rjr] and KK final states as well. This paper presents information which may be used in 
such a program. 

The J^'-^ of the 7r'^7r'^ system must have J even with both P and C positive. The isospin must also be 
even (I = or I = 2) for 7r''7r*^. The 7r*^7r*^ system is thus particularly attractive for investigation of scalar 
and tensor states as the PWA is simplified without the presence of odd angular momenta. 

The Tr~p — > Tr^w^n reaction has been studied in experiments with incident 7r~ momenta of 9 GeV/c [ p^ , 
25GeV/c 38GeV/c [|l9| and 100 GeV/c The combined information from these experiments can 

be used to provide information on how cross sections of produced states and relative ratios of partial waves 
depend on center-of-mass energy. 

This paper is organized as follows: The experimental overview is presented in Section 2. Event recon- 
struction and data selection are described in Section 3, where the general features of the distributions in 
rn^o^o and \t\ are also discussed. The details of the PWA formalism and results are given in Section 4. In 
Section 5 Regge-models are fitted to the results from Section 4. The conclusions are summarized in Section 
6. 



II. EXPERIMENTAL OVERVIEW 

The E852 apparatus |2^ was built around and included the Multi-Particle Spectrometer (MPS) at 
BNL. The data used for the analysis reported in this paper were collected in 1994 and 1995 using a beam 
of negatively charged particles of momentum 18.3 GeV/c. A 30-cm liquid hydrogen target was surrounded 
by a cylindrical drift chamber |Q and an array of thallium-doped Csl crystals arranged in a barrel, 
all located inside the MPS dipole magnet. Drift chambers were used to track charged particles downstream 
of the target. Two proportional wire chambers (PWC's), downstream of the target, were used in requiring 
specific charged particle multiplicities in the event trigger. A 3000-element lead glass detector (LGD) ||2^ 



2 



measured the energies and positions of photons in the forward direction. The dimensions of the LGD matched 
the downstream aperture of the MPS magnet. Photons missing the LGD were detected by the Csl array or 
by a lead/scintillator sandwich array (DEA) arranged in a picture frame downstream of the target with an 
aperture to allow for the passage of charged particles. 

The first level trigger required that the unscattered or elastically scattered beam not enter an arrange- 
ment of two small beam-veto scintillation counters located in front of the LGD. The next level of trigger 
required that there be no signal in the DEA and no charged particles recorded in the cylindrical drift chamber 
surrounding the target or in the PWC's (an all- neutral trigger). In the 1994 run, all layers of the cylindrical 
drift chamber were used in the trigger requirement, whereas in the 1995 run, only the outer layer was used. 
A common off-line analysis criterion required no hits in the cylindrical drift chamber. The final trigger 
requirement was a minimum deposition of electromagnetic energy in the LGD. 

The LGD is central to this analysis and it is described in detail in reference |Q. The LGD was initially 
calibrated by moving each module into a monoenergetic electron beam. Further calibration was performed 
by adjusting the calibration constant for each module until the width of the 7r° and ry peaks in the 77 effective 
mass distribution was minimized. The calibration constants were also used for a trigger processor which did 
a digital calculation of energy deposited in the LGD and the effective mass of photons striking the LGD. A 
laser-based monitoring system allowed for tracking the gains of individual modules. 

Studies were made of various algorithms for finding cluster of energies deposited by photons including 
issues of photon-to-photon separation and position finding resolution. These are also described in reference 




III. EVENT RECONSTRUCTION AND DATA SELECTION 



The combined data sets taken in 1994 and 1995 contain approximately 70 million all- neutral triggered 
events. Of these events, approximately 13 million were found to have four photons in the LGD. The plot 
of di-photon effective masses for all possible pairings of photons is shown in the scatterplot of figure |l|(a) 
and the projection is shown in figure |l|(b). Events consistent with the production of two tt^'s dominate the 
scatterplot. The tt" mass resolution is 17 MeV/c^ . The sample of 847,460 ■K~p Tr^Tr^n events was selected 
from the 13 million four photon events by imposing various analysis criteria. It was required that no charged 
particles were registered in the MPS drift chambers or the cylindrical drift chamber surrounding the liquid 
hydrogen target. Any event with a photon within 8 cm of the center of the beam hole or the outer edge of 
the LGD was removed. The returned from kinematic fitting to the np —> Tr'^Tr'^n reaction hypothesis was 
required to be less than 9.8 (95% G.L. for a three-constraint fit). A further demand was that none of the 
other final state hypotheses considered {ijir'^n, ijrjn) had a better x^- The final criterion was that the Gsl 
detector registered less than 20 MeV, a cut which eliminated events with low-energy tt'^'s. The 7r°7r° mass 
resolution improves from 24 MeV/c^ to 16 MeV/c^ at the mass of the after kinematic fitting. 

Background studies were also carried out. By selecting events in a given four photon effective mass region 
and fitting the associated scatter plot of di-photon effective mass pairings (similar to figure ||), the background 
of non-Tr'^Tr'' events under the signal was found to be very small. Typical signal to noise ratios determined 
by these studies are in the range of 50:1. Monte Carlo studies indicate that combinatoric background from 
mis-pairing the reconstructed photons is a few percent below m^jr ~ 0.5 GeV/c^ and non-existent at higher 
masses. These studies are described in more detail in reference |p5| . 

The distribution in missing-mass-squared, recoiling against the four photons, for events with a successful 
kinematic fit to the reaction n^p — > Tr^Tr^n is shown in figure The missing-mass-squared is determined 
from photon position and energy information before kinematic fitting and the distribution peaks near the 
square of the neutron mass. The distribution in tt^tt^ effective mass is shown in figure ^ The spectrum is 
dominated by the /2(1270) resonance and a broad enhancement at low tt'^tt^ mass (from threshold to about 
1.0 GeV/c^). There is also a small Kg tt^tt^ signal present, despite the requirement that the deposited 
energy in the Gsl detector not exceed 20 MeV. This Gsl energy cut reduces a substantial fraction of Kg 
events but other reactions producing Kg can avoid deposition of energy in the Gsl detector. By correlating 
the observed yields of Kg and /2(1270) mesons, for samples with and without the Gsl detector energy cuts, 
with cross sections for /2(1270) production and associated Kg production {n^p KgA{I]'^)) measured in 



3 



other experiments, we estimate an overall Csl detector inefficiency of 5%. These studies also indicate that 
the background level of non-neutron events under the /2(1270) is approximately 1%. Another feature of the 
spectrum is the dip at 1.0 GeV/c^ , which will be seen to be due to the interference of a narrow resonance, 
the /o(980), with a broad tt'^tt^ enhancement. 

The distribution in |t|, shown in figure |[ is not characterized by a single exponential, suggesting more 
than one production mechanism. The curve is a fit of this distribution to a sum of two exponentials: 
dN/dt — a-e~'' l*l +c-e^'* l*l where b = 15.5 {GeV/c)~'^ and d = 3.7 (GeF/c)~^. Based on this structure, we 
initially examine the tt^tt^ effective mass spectra in four bins in \t\ as shown in figure |[ The i-dependence 
of the S, Dq, and D+ partial waves is later investigated in a set of partial wave fits more finely binned in \t\. 

An inspection of figure ^ reveals striking differences in the tt^tt^ mass spectra associated with the four 
bins in \t\. For example, the low-mass structure which dominates in figure ^ is much less prominent in 
figure The dip associated with the /o(980) resonance in figure ^ becomes a bump in figure These 
and other features are explored in more detail below in the discussion of the PWA results. 



IV. PARTIAL WAVE ANALYSIS 



Partial wave analysis is used to extract production amplitudes (partial waves) from the observed decay 
angular distributions of the di-pion system. A process such as n^p — > Tr'^Tr'^n, dominated by t-channel meson 
exchange, is simplest to analyze in the Gottfried- Jackson reference frame. The Gottfried- Jackson frame is 
defined as a right-handed coordinate system in the center of mass of the produced di-pion system with the 
z — axis defined by the beam particle momentum and the y — axis perpendicular to the plane defined by 
the beam and recoil neutron momenta. The decay angles {6, (p) are determined for one of the produced 7r° 
momenta. At fixed beam momentum, an event is fully specified by {mT^Tr,t,9, (j)). The data are binned in 
m^TT and t and the production amplitudes, and their relative phases, are extracted from the accumulated 
angular distributions using an extended maximum likelihood fit to the distributions in (0,0) [ p6| . The naming 
convention for the partial waves is summarized in Tabic |[ 

The explicit form of the angular distribution I{9, (j)) fitted to the data in a given mass and momentum 
transfer range in this analysis is given by 



I{9. 



S + V5DoP§{cost 



cos 9) cos i; 



- V9GoPi'(cos6') 



D+P^ {cos 9) sin( 



(1) 



where P["{9) are the associated Legendre polynomials |2q ]. 

As summarized in Table the D+ wave is produced by the exchange of a particle with natural parity 
(P = (— l)''). For production of a rnr system, the dominant natural parity exchange particle is the 02 
p7| . The S, Do, D- and Go waves are produced by the exchange of a particle with unnatural parity 
(P — (— l)"'^^). Again, for tttt production, the dominant unnatural parity exchange particles are the tt and 
the Oi p7|. 



A. Ambiguities 

There are multiple discrete sets of partial wave amplitudes which can give rise to exactly the same 
angular distribution ]2^ ]. It can be shown that in a partial wave fit with only S, Dq, D^, and 13+ partial 
waves there are four sets of ambiguous partial wave amplitudes. The four sets can be divided into two 
groups with different partial wave intensities. Additionally, within each group, there is a sign ambiguity in 
the phases between the amplitudes. 

Normally there are two ambiguities; if one wave in each naturality is fixed a priori, e.g. set real or to 
some complex value for dynamical reasons, there is still an overall sign ambiguity. However, even this sign 
ambiguity could be fixed by the requirement, for example, of a resonant behavior in one of the waves. 
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In general, there is an eightfold ambiguity for a tt^tt^ system containing L=0, 2 and 4. However, these 
ambiguities necessarily entail nonzero G_ and G+ waves. In this paper we have assumed that these are 
negligibly small and searched for ambiguities with nonzero Go wave. We find no such ambiguities in our 
data. 

In the analysis of the tt^tt*' system, the physical solution can be selected by a combination of physical 
arguments (which will be given below) and the requirement that solutions be smoothly connected as a 
function of mass. This selection of the physical solution applies simultaneously to all intensities and phases. 
In what follows, the physical solution is plotted with solid symbols. The other solutions are plotted with 
open symbols and are presented for completeness. 



B. Partial Wave Fits 



1. Results for 0.01 < -t < 0.10 GeV'^/c^ 

The results of the partial wave decomposition are shown in figures ^ and ^. The partial wave intensities 
are shown in figure]^ and the phase differences in figure^. The phase difference plots are shown above tt^tt^ 
masses of 0.8 GeV/c^. Below that value, where one of the waves is very small, phase difference information 



is unreliable. As discussed in Section IVA. there is a two- fold ambiguity in the intensities. The threshold 
behavior (S'-wave dominance) and the resonant behavior of the /2(1270) are used to select the physical 
solution. Furthermore, since the resonant structures of both the Dq and -D_ partial waves are due to the 
/2(1270), the relative phase between the Dq and D_ partial waves should be constant and near ±7r radians, 
according to the phase convention of [^6| . These assumptions allow the physical solution at low mass to 
be connected with the solutions at higher mass. Above approximately 1.5GeV/c^, the solutions become 
degenerate. The spin-4 Gq partial wave is not included in the fit below lAGeV/c^ . 

There are a number of key features observed in the physical solution. There is at least one broad 
enhancement in the S'-wave intensity and a sharp dip in the S'-wave intensity near 1.0 GeV/c^ accompanied 
by rapid phase variation in the S— Dq relative phase. There also exists a dip in the S'-wave intensity near 
1.5 GeV/c^ accompanied by rapid phase variation in the S— Dq relative phase. The /2(1270) is observed in 
the Dqj and partial wave intensities, and the bump observed in the Go partial wave near 2.0 GeV j (? 
is consistent with the /4(2040). Finally, the Do-wave intensity is larger than the I?_-wave intensity or the 
ZJ^-wave intensity, consistent with the expectation that OPE should favor production of an m=0 wave for 
this low-|f| region. 

A background term was not included in the PWA fits presented in this paper. A background term was 
included in some earlier fits where it was found that below about 1.0 GeV j (? it cannot be distinguished from 
the dominant S'o wave and above 1.0 GeVjf? , the fit forces the background term to zero. 



2. Results for 0.10 < -t < 0.20 GeV^/c^ 



The results of the partial wave analysis (Figures ^ and ^ in this region are qualitatively similar to the 
results in the 0.01 < —t < 0.10 GeV^ /c^ region. The same techniques are used to select the physical solution 
as in the previous region in The S'-wave intensity contains at least one broad object and two dips. The 
/2(1270) is observed in all D-waves. An enhancement near 2.0 GeV/c^ is again observed in the Gq partial 
wave. More detailed comparisons with the results from the 0.01 < — t < 0.10 GeV/c^ region reveal the 
following differences: The ratio of the S'-wave intensity to the I?o-wave intensity is smaller at larger \t\ and 
the ratio of the Z?o-wave intensity to both the Z?_-wave and Z?+-wave intensities is smaller at larger \t\. The 
ratio of the Do-wave intensity to Go-wave intensity does not change suggesting that the /2(1270) and the 
/4(2040) are produced by the same mechanism. 
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3. Results for 0.20 < -t < OAOGeV'^/c^ 



The change in slope for the |i|-distribution as seen m figure indicates a change in production mech- 
anism. This is reflected in the partial wave analysis as well (Figures ^ and ^0|). For this \t\ region, the Go 
partial wave is not required for an adequate description of the observed angular distributions and is therefore 
not included here or in the next higher \t\ region. The S'-wave intensity has a different shape compared to 
that at smaller values of \t\. The Z?+-wave intensity is approximately one-third as large as the Z^o-wave 
intensity whereas at smaller momentum transfer it was approximately one-tenth as large. 

4. Results for 0.40 < -t < 1.50 GeVVc^ 

The partial wave analysis results in the region 0.40 < —t < l.bGGeV^ / (Figures |ll| an djl^ ) are 
significantly different from results at smaller The bump observed in the mass plot (figure ^d) near 
1.0 GeVji? is found in the 5*- wave intensity. The -D+ partial wave is dominant (as opposed to the Dq partial 
wave), indicating a shift from unnatural parity exchange processes at small |t| to production via natural 
parity exchange at large \t\. 

5. Fine \t\ Bin Fits 

The statistics of this experiment are sufficient to allow the region 0.00 < —t < OAOGeV^/c^ to be 
analyzed in finer \t\ bins, nine in all, for masses up to approximately 1.8 GeV/c^. The \t\ -dependence 
of the S'-wave intensity may be summarized by noting that the ratio of the maxima in the intensities at 
approximately 0.8 GeV/c^ and 1.3GeV/c^ decreases with increasing \t\, and the ratio of height of maximum 
intensity at approximately 0.8 GeV/c^ to value of the intensity measured at 0.98 GeV/c^ decreases. 

The lineshape of the /2(1270) in the Do-wave intensity is largely independent of \t\. The S — Dq relative 
phase is |t| dependent. The lineshape of the D_|.-wave is also independent of |t|. More details of the |i| 
dependence of the partial waves follow. 

The intensities of the individual partial-waves and phase differences as a function of mass for the nine 
bins in |t| for 0.00 < —t < 0.40 GeV'^/c^ as well as for the |t|-bins presented in this paper are available on 
the World Wide Web !§. 

C. Model Dependent Fits of the \t\- Distributions 

The integrals of fitted relativistic Breit-Wigner functions over the peak regions of the Dq and D^-waves 
as a function of |i| are shown in Figure The dependences of these intensities on |i| are fitted to functions 
given by Regge-exchange models. At low-|i|, the unnatural parity exchange Dq partial wave is expected to 
be dominated by OPE. The Reggeized form for this contribution is given by 

^ = iVz,jy=te^'*(t-m2j2(l-f e-"(*))r(-a.(i))P, a.{t) ^ 0.9{t - ml) (2) 

In this expression, the factor is due to helicity-flip in the pion-nucleon coupling, and the polynomial 
dependence on t arises from the /2 coupling to tttt at the production vertex. The particular form of this 
dependence is due to the angular momentum barrier factor proportional to fc^ with L — 2 and k being the 
magnitude of the 3-momentum of the exchanged particle in the /2 rest frame (Gottfried- Jackson frame), 
given by fc^ = ((to/, - to^)^ - t){{mf^ +m^f - t)/Amj^ ~ {mj^ - tf / 'Im)^. 

The slope, 6^, in the OPE form is 4.08±0.02/(Gey^/c^). The systematic uncertainty in the slope of the 
a^(i) Regge trajectory is ±0.1/(Gey^/c^). As shown by Irving and Michael |§] the natural parity exchange 

-wave is dominated by absorption of the pion exchange and may be parameterized in terms of a Regge 
cut in the nucleon helicity-flip amplitude 
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ac(t)-l 

C = ffee'-^'e-^^^^^^*) ^ , ac{t) = 0.41i, = -0 .84, 6, = 3.89 (3) 



The nucleon-flip and non-flip 02 exchange is then given by, 



respectively, with the parameters a^glO ~ 0-5 + 0.82t and = 1.35, ha — 3.24, po — 17.2 Gey/c, and 
r = 0.5 from jsj and — 18.3 GeV/c^ the beam momentum for these data. The _D+-wave intensity is then 
fitted to 

d\D+ 



d\t\ 



ND^{\An\^ + \Af + C\^) (5) 



For both forms the fitted functions are averaged over the \t\ bins shown in the plots. The plotted curves 
are calculated from the models without averaging. 

In figure ImI the peak value of the S-wave intensity near 0.80 GeV/c^, the value of the S'-wave intensity 
at 0.98 Ge^/c^^ and the peak value of the S'-wave intensity at approximately l.iGeV / as a function of 
\t\ are shown. A one-pion-exchange form similar to Equation ^, but with the t — rri^f^ factor removed, is 
overlayed on these distributions. The Regge trajectory slope and exponential slope are fixed to the values 
found for the Dq-wsnb fit, and a one parameter fit is used to set the normalization. At small values of 
\t\ the OPE form qualitatively agrees with the data. The excess of events at higher |t| in (b) and (c) is 
consistent with the existence of additional production mechanisms that are less strongly biased toward small 
momentum-transfer-squared production than is OPE. 



V. CONCLUSIONS 



A partial wave analysis was carried out on a sample of 847,460 events of the reaction n^p ir^TT^n 
collected by experiment E852. The PWA was performed in 0.04GeV^/c^-wide bins in di-pion mass (to^o^o) 
and momentum-transfer-squared (|<|) from the incident tt^ to the outgoing tt^tt*^ system. Coarse and fine 
binning in |i| were used. Numerical values for the partial wave intensities and phases as a function of di-pion 
mass, for coarse and fine bins in |t|, are available on the World Wide Web The /2(1270) meson is found 
to be produced by unnatural parity exchange at small values of |f| and natural parity exchange at large 
values of \t\. The |t| dependences of I?o-wave and D+-wave intensities are consistent with Regge-exchange 
models. An enhancement in the Gq wave consistent with the /4(2050) meson is observed in unnatural parity 
exchange at small momentum transfer. The shape of the 5-wave intensity has a strong momentum transfer 
dependence. The presence of dips in the S'-wave intensity near 0.98 and 1.5 GeTZ/c^, accompanied by rapid 
phase variations relative to the Dg-wave is consistent with similar observations reported in reference ||2(][| 
and in centrally produced tt^tt^ systems in 450 GeV/c pp collisions ||2^. The latter claims evidence for the 
/o(980) and /o(1500). At large momentum transfer, the /o(980) meson is observed as a bump in the S-wave 
intensity. The S-wave intensity in the peak near 0.80 GeV/c^ is well-described by OPE. It should also be 
noted that the model of Anisovich et al ||ll|-[l3|| predicts the presence of a dip in the \t\ distribution for this 
mass region near \t\ « 0.07Gey^/c^ [Q. In direct contradiction, no such dip is observed in this analysis. 
At higher masses the S-wave is adequately described by OPE only at small values oi\t\. 
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FIG. 1. (a) The plot of pairs of di-photon effective masses {rriij vs. rriki) for all pairs of photons {i,j,k,l) is 

dominated by the Tr^n° signal. Clear evidence is also seen for the production of rj-n-". (b) The projection of the scatter 
plot is shown. 
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FIG. 2. The missing-mass-squared distribution is shown for four-photon events with a successful kinematic fit to 
the reaction ■K~p —^ 7r^7T°n 
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FIG. 3. (a) The n^n'^ effective rticiss distribution. This spectrum is dominated by the presence of the /2(1270) 
resonance. Additionally, there is a broad enchancement peaking near 0.8 GeV/c^ and dips in the spectrum at 1.0 

and 1.5 GeV/c^. (b) The momontum-transfor-squarcd distribution with a fit to the sum of two exponential functions. 
The structure of this distribution is suggestive of changing production mechanisms. 
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FIG. 4. The tt'V" effective mass distribution for four regions of \t\. The shape is seen to be strongly dependent 
on \t\. Of particular interest is the disappearance of the broad enhancement near 0.8 GeV/c^ in (a) and (b) and the 
emergence of a small peak at 0.98 GeV/c^ in (c) and (d) with increasing values of \t\. 
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FIG. 5. The squares of the magnitudes of the partial waves (a)-(e) as a function of mass for events in the region 
0.01 < \t\ < Q.lQGeV^ /(? . The solid circles correspond to the physical solution. The open circles correspond to 

the other ambiguous solution. Additionally the coherent sum of the partial waves integrated over decay angles, (f), 
gives the acceptance corrected mass distribution. The dominant production mechanism is m = unnatural parity 
exchange {S , Do , and Go partial waves). 
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FIG. 6. For events in the region 0.01 < \t\ < 0.10 the relative phase between the S and Dq partial waves (a) shows 
rapid phase variation near 0.98 GeF/c^ and 1.5 GeV/c^. The relative pi lase between the Do and D_ partial waves 
(b) is smooth and nearly constant up to 1.5 GeV/c^. The relative phase between the S and Go partial waves is shown 
in (c). The solid circles represent the physical solution. 
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FIG. 7. The squares of the magnitudes of the partial waves (a)-(e) as a function of mass for events in the region 
0.10 < \t\ < 0.20 GeV^ /c^ . The sohd circles correspond to the physical solution. The coherent sum of the partial 
waves integrated over decay angles, (f), gives the acceptance corrected mass distribution. As in figure ^, the dominant 
production mechanism is m = unnatural parity exchange {S , Do , and Go partial waves). 
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FIG. 8. For events in the region 0.10 < \t\ < 0.20 the relative phase between the S and Do partial waves (a) shows 
rapid phase variation near 0.98 GeF/c^ and 1.5 GeV/c^. The relative pi lase between the Do and D_ partial waves 
(b) is smooth and nearly constant up to 1.5 GeV/c^. The relative phase between the S and Go partial waves is shown 
in (c). The solid circles represent the physical solution. 
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FIG. 9. The squares of the magnitudes of the partial waves (a)-(d) as a function of mass for events in the region 
0.20 < |t| < QAQ GeV^ / (? . The solid circles correspond to the physical solution. The coherent sum of the partial 
waves integrated over decay angles, (e), gives the acceptance corrected mass distribution. Compared to figures |^ 
and ^ the D+ partial wave (natural parity exchange) is becoming more important although the dominant production 
mechanism is still m = unnatural parity exchange [S , Do , and Go partial waves). 
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FIG. 10. The relative phases between unnatural parity exchange partial waves for events in the region 
0.20 < \t\ < 0.40 . The physical solution (solid circles) in the S — Do relative phase plot (a) shows less rapid 
phase variation than in figures ^ and|^. The Do — D- relative phase is shown in (b). 
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FIG. 11. The squares of the magnitudes of the partial waves (a)-(d) as a function of mass for events in the region 
0.40 < \t\ < 1.50 GeV^ /c^ . The solid circles correspond to the physical solution. The coherent sum of the partial 
waves integrated over decay angles, (e), gives the acceptance corrected mass distribution. The D+ partial wave 
(natural parity exchange) is the dominant partial wave. 
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FIG. 12. The relative phases between unnatural parity exchange partial waves for events in the region 
0.40 < \t\ < 1.5 . The physical solution (solid circles) in the S — Do relative phase plot (a) is more smoothly 
varying than in figures ^ and ^. The Do — D- relative phase (b) is constant only up to approximately 1.2 GeV/(? . 
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FIG. 14. The S-wave intensities at three different masses (0.80, 0.98, and 1.30 GeV/(?) are compared with a 
one-pion exchange form. Except for the overall normalization, the parameters of the OPE parameterization are those 
determined in the fit to figure ^(a). OPE describes the data well at small \t\. The excess of events at higher \t\ in 
(b) and (c) is consistent with the existence of additional production mechanisms that are less strongly biased toward 
small momentum-transfer-squared production than is OPE. 
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Partial wave 


L 


\m\ 


Naturality of the exchange particle 


S 








unnatural 


Do 


2 





unnatural 


D- 


2 


1 


unnatural 


Go 


4 





unnatural 


D+ 


2 


1 


natural 



TABLE I. The nomenclature for partial waves includes the angular momentum (L) of the 7r°7r° system, the 
magnitude of the magnetic quantum number (m) and the naturality of the exchange particle which leads to production 
in the particular partial wave. The naturality is natural if P = (—1)'' and unnatural if P = (—1)'^"'"^. 
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